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Abstract. We present relations between X-ray luminosity and velocity dispersion (L - a), X-ray luminosity and gas mass 
(L - A/ ga8 ), and cluster radius and velocity dispersion (r 50 o - c) for 62 galaxy clusters in the HIFLUGCS, an X-ray flux-limited 
sample minimizing bias toward any cluster morphology. Our analysis in total is based on ~ 1.3 Ms of clean X-ray XMM-Newton 
, data and 13439 cluster member galaxies with redshifts. Cool cores are among the major contributors to the scatter in the L- cr 

relation. When the cool-core-corrected X-ray luminosity is used the intrinsic scatter decreases to 0.27 dex. Even after the X-ray 
luminosity is corrected for the cool core, the scatter caused by the presence of cool cores dominates for the low-mass systems. 
The scatter caused by the non-cool-core clusters does not strongly depend on the mass range, and becomes dominant in the 
high-mass regime. The observed L-cr relation agrees with the self-similar prediction, matches that of a simulated sample with 
AGN feedback disregarding six clusters with < 45 cluster members with spectroscopic redshifts, and shows a common trend of 
increasing scatter toward the low-mass end, i.e., systems with <x < 500km s -1 . A comparison of observations with simulations 
indicates an AGN-feedback-driven impact in the low-mass regime. The best fits to the L - M gas relations for the disturbed 
clusters and undisturbed clusters in the observational sample closely match those of the simulated samples with and without 
AGN feedback, respectively. This suggests that one main cause of the scatter is AGN activity providing feedback in different 
phases, e.g., during a feedback cycle. The slope and scatter in the observed r 50 o - cr relation is similar to that of the simulated 
sample with AGN feedback except for a small offset but still within the scatter. 
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d ' 1 . Introduction During the past, there have been a large number of studies of 

X-ray luminosity scaling relations along with their applications 

Galaxy clusters have been suggested as a potential probe of the tQ confjtrain cosmological parameters in ga i axy c i uster surveys 

dark energy equation of state parameter (w = p/p, where p is and ±s physical ^ of ±e hot intracluster medium (ICM) in 

the energy density and p is the pressure), through the evolu- galaxy dusters (g g ^ Henry & Tucker mg . Hemy & Amaud 

tion of the mass function (e.g., Schuecker et al. 2003; Predehl mi . Edge & Smt ml . Dayid e( a , m3 _ Fabian gt aJ 

et al. 2007; Henry et al. 2009; Vikhlinin et al. 2009a, 2009b; m4 . Girardi e( a , m& Mushotzky & Scharf 1997; Ca valiere 

Mantz et al. 2010). Observational surveys select galaxy clus- gt ^ mq . whhe gt aJ m7 . Markevitch 1998; Wu et al 

ters by their observables rather than by their mass. Therefore, 199g> 1999; ^ & Fabian 199g . Amaud & Eward 1999; 

a relationship between the cluster total mass and an observ- Reiprich & B5hringer 2002; Qta et al 2006; chen et al 2007; 

able such as X-ray luminosity is required to recover the se- zhang rf aJ 2006i 2008; praU gt aJ 2QQ9 . Leauthaud et al 

lection function of an X-ray survey in terms of cluster masses 201Q . Stanek ^ al mQy Lafge x _ ray cosmology surveySj e g ^ 

and predict the cluster mass, hence the cluster mass function. fey eR0SITA m expected to substantially improve cosmolog- 
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eROSITA, the use of X-ray mass proxies has been proposed, 
specifically X-ray luminosity, to infer the total mass and con- 
struct the selection function in the forthcoming wide survey of 
the satellite (Predehl et al. 2007). The superb quality X-ray data 
in the XMM-Newton archive provide us with an excellent op- 
portunity to calibrate the luminosity scaling relations and more 
clearly understand the X-ray selection method. 

Simulations show that the formation of galaxy clusters is 
not a purely gravitational process; The galaxy velocity dis- 
persions of clusters appear to indicate that heating is present 
when compared to the cold dark matter (CDM) velocity disper- 
sion normalized to the Wilkinson Microwave Anisotropy Probe 
(WMAP) and large-scale structure (LSS) distributions (Evrard 
et al. 2008). Cluster mergers not only change the cluster X-ray 
luminosity (e.g., Ricker & Sarazin 2001; Poole et al. 2006), but 
also affect the properties of the cluster galaxies (e.g., Sun et al. 
2007; Smith et al. 2010). Although the hot gas and galaxies are 
not pure tracers of the gravitational potential of galaxy clusters, 
they are indeed sensitive probes of the dynamical properties 
of galaxy clusters, and react on different timescales during a 
merger in simulations (e.g., Roettiger et al. 1999). In particular, 
the optical information about the line-of-sight velocity of clus- 
ter galaxies complements X-ray information about the cluster 
morphology projected onto the sky. The luminosity - velocity 
dispersion (L - cr) relation of galaxy clusters is thus crucial to 
understanding the dynamical properties of galaxy clusters and 
their impact on the scaling relations and possibly the X-ray se- 
lection bias (e.g., Wu et al. 1999; Ortiz-Gil et al. 2004). 

To carry out the L-cr studies, one requires a representative 
sample with a well-defined selection function and minimal bias 
toward any cluster morphology, as well as superb quality X-ray 
data and large amount of cluster galaxy redshifts. The Highest 
X-ray FLUx Galaxy Cluster Sample (HIFLUGCS, Reiprich 
& Bohringer 2002) of 64 galaxy clusters selected from the 
ROSAT All-Sky Survey (RASS; Ebeling et al. 2000; Bohringer 
et al. 2004) is such a sample. In the HIFLUGCS, we analyzed 
all available X-ray data in the XMM-Newton archive for 63 
clusters which represents nearly 4 Ms of data. After clean- 
ing and selecting the longest observation closest to the clus- 
ter center for clusters with multiple observations, we still have 
~ 1.3 Ms XMM-Newton data for 59 clusters. For 62 clusters 
in the HIFLUGCS, we obtained a sum of 13439 cluster mem- 
ber galaxies based on spectroscopic redshifts and performed a 
careful exclusion of non-members. In the end, we were able 
to measure X-ray observables, combining XMM-Newton and 
ROSAT data, and velocity dispersion, based on 13439 cluster 
members, to make a cross-calibration for 62 out of 64 clusters 
in the HIFLUGCS. 

The outline of this paper is as follows. We describe the data 
analysis in Sect. [2] present the scaling relations of the 62 clus- 
ters in the HIFLUGCS in Sect. [3] compare the observational 
and simulated samples in Sect. |4] discuss the systematic er- 
rors in determining the velocity dispersion in Sect. [5] and sum- 
marize our conclusions in Sect. [6] Our Appendix provides ex- 
tra information on the cross-calibration between XMM-Newton 
and ROSAT, the iron abundance versus (vs.) temperature cor- 
relation, results using either the 0.5-2 keV X-ray luminosity 
corrected for the presence of a cool core (< 0.2r 5 oo), or the 



luminosity including or excluding the cluster core, the XMM- 
Newton 0.7-2 keV images, and the figures illustrating system- 
atic errors in estimates of cr. Throughout the paper, we as- 
sume that Q m = 0.3, Q A = 0.7, and H = 70 kms -1 Mpc -1 . 
Confidence intervals correspond to the 68% confidence level. 
Unless explicitly stated otherwise, we apply the BCES regres- 
sion fitting method taking into account measurement errors in 
both variables (Akritas & Bershady 1996). 

2. Data analysis 

2.1. Optical data analysis and velocity dispersion 

We draw the velocity of the cluster galaxies from the litera- 
ture (updated until March 2010, including the compilation in 
Andernach et al. 2005). When there is more than one veloc- 
ity per galaxy, we calculate an averag^] of the measurements, 
excluding discordant values and those with large errors when 
more than one measurement is available. 

Brightest cluster galaxies (BCGs) in galaxy clusters are al- 
most invariably giant ellipticals and are more luminous than 
normal galaxies. The BCGs have line-of-sight velocities that 
are similar to the mean of their host clusters and extended 
stellar envelopes. We identify the BCG on the basis of its ap- 
parent magnitude and spectroscopic confirmation as a cluster 
member. To define a BCG position for every HIFLUGCS clus- 
ter, we made the following choices for clusters without a sin- 
gle dominant BCG. A0400 and A2065 have dumbbell BCGs, 
and A3 158 and A2256 have BCG pairs, for which we place 
the BCG positions in the middle of the two components of 
indistinguishable brightness. A3266, A3391, A0576, A2634, 
MKW8, and IIIZw54 have dumbbell BCGs, and Coma and 
Hydra (A1060) have two brightest galaxies of similar bright- 
ness, for which we place the BCG position on the brighter 
component as the difference in the brightness is measurable. 
A2199 has multiple nuclei, and we place the BCG position at 
the brightest nucleus. We list the BCG positions in TableQ] 

As most BCGs are located very near the X-ray flux- 
weighted cluster centers (definition see § 12. 3. Il l, we select pre- 
liminarily galaxies with spectroscopic redshifts in each cluster 
within an aperture of at least 1.2 Abell radius, i.e., 2.57 Mpc, 
centered on the BCG. For each cluster, we plot the line-of- 
sight velocity of the selected galaxies as a function of their 
projected distance from the BCG, and locate the caustic, a 
trumpet-shaped region, which efficiently excludes interlopers 
(e.g., Diaferio 1999; Katgert et al. 2004; Popesso et al. 2005; 
Rines & Diaferio 2006). We consider only the galaxies inside 
the caustic as cluster members, and exclude the others from 
subsequent analysis. More than 80% of the clusters have a 
clearly evident caustic shape. In Fig. [T] we show as an ex- 
ample the caustic and sky positions of the cluster galaxies in 
Coma, and as a poor example in SI 101. There are eight clus- 
ters that have fewer than 45 cluster members with spectro- 
scopic redshifts in the HIFLUGCS. We excluded 2A0355 and 
RXCJ1504 from our study since both have at most three red- 
shifts each. We still consider the remaining six systems (i.e., 

1 Since the individual error estimates are inhomogeneous, we de- 
cided not to weight the calculation of the average velocity. 
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A0478, NGC1550, EXO0422, HydraA, SI 101, and A2597) 
with > 12 but < 45 cluster members with spectroscopic red- 
shifts, and highlight them in our results. We gathered a total 
of 13439 cluster- member galaxies based on spectroscopic red- 
shifts and a careful exclusion of non-members, which gives a 
median value of 185.5 per cluster. 

For the galaxies selected as cluster members we apply the 
bi-weight estimator (e.g., Beers et al. 1990) to measure the ve- 
locity dispersion. The errors are estimated through 1000 boot- 
strap simulations. We list the number of cluster members (« ga i) 
and the velocity dispersion (<x) of the cluster for 62 clusters in 
the HIFLUGCS in Table Q] The systematic errors in the deter- 
mination of the velocity dispersion are discussed in Sect. [5] 

2.2. X-ray data analysis 

There are 63 clusters in the HIFLUGCS in the XMM-Newton 
archive. Only A2244 has not yet been observed. We ana- 
lyzed 150 XMM-Newton observations, which give 3.90 Ms for 
MOS1, 3.97 Ms for MOS2, and 3.68 Ms for pn, respectively. 
To filter flares, we apply iterative screening similar to Zhang et 
al. (2006) using both the soft (0.3-10 keV) band and the hard 
(10-12 keV for MOS, 12-14 keV for pn) band but with a 3.3- 
cr clipping. We found that the XMM-Newton observations of 
four clusters (i.e., A0401, A0478, A1736, A2163) are flared. 
For clusters with multiple observations, we select the longest 
observation of which the pointing position is the closest to the 
cluster center. Since 2A0355 and RXCJ1504 have at most three 
redshifts each, we exclude these two clusters, and end up with 
nearly 1.3 Ms of XMM-Newton data of 57 clusters for a more 
detailed analysis. 

For the four flared clusters (i.e., A0401, A0478, A1736, 
A2163), as well as for A2244, the X-ray quantities are de- 
rived from ROSAT pointed observations. The X-ray quantities 
for the remaining 57 clusters in the HIFLUGCS are derived 
from combined XMM-Newton and ROSAT data. We note that 
the XMM-Newton observations only cover an incomplete sec- 
tor of A2142, such that we have to use the ROSAT data to de- 
rive its surface brightness profile. The XMM-Newton observa- 
tions of A2142 are only used to measure the global temperature 
and iron abundance. We describe in detail the procedures we 
adopted to detect and subtract point-like sources and for back- 
ground treatment in Sect. 2 and Sect. 3 of Zhang et al. (2009). 
Significant substructure features clearly detected in the image 
are excised before we perform the spectral and surface bright- 
ness analysis. We note that the surface brightness analysis is 
slightly different from that in Sect. 4 of Zhang et al. (2009) in 
that we directly convert the ROSAT surface brightness profile 
to the XMM-Newton count rate using the best-fit spectral model 
obtained from the XMM-Newton data. We then combine the 
XMM-Newton surface brightness profile within the truncation 
radius, where the XMM-Newton signal-to-noise ratio is ~ 3, 
with the ROSAT converted surface brightness profile beyond 
the truncation radius for further analysis. We list the properties 
of the XMM-Newton observations, redshift, hydrogen column 
density, gas mass, X-ray morphology, and presence of a cool 
core of each cluster in Table [2] 



2.2.1 . Cluster radius, i.e., r 50 o 

X-ray quantities have to be derived consistently within a certain 
cluster radius, e.g., rsoo, the radius within which the mass den- 
sity is 500 times of the critical densit}0, at the cluster redshift. 
The quantity rsoo can be measured from the X-ray measured 
mass distribution derived under the assumption of hydrostatic 
equilibrium as we did in Zhang et al. (2009). Observations 
have found evidence of deviations from hydrostatic equilib- 
rium (e.g., Zhang et al. 2008, 2010; Mahdavi et al. 2008). The 
cross-calibration between weak lensing masses and X-ray ob- 
servables instead uncovers a tight scaling relation between gas 
mass and cluster total mass (e.g., Okabe et al. 2010). We there- 
fore use the gas mass to infer the cluster mass and rsoo- Our 
sample occupies a wide mass range with the gas masses from 
1.74 x 10 n M o to 2.12 x 1O 14 M , which is similar to the mass 
range of the extended sample in Pratt et al. (2009) consisting of 
41 groups and clusters collected from Vikhlinin et al. (2006), 
Arnaud et al. (2007), Bohringer et al. (2007), and Sun et al. 
(2008). We thus adopt their relation E h5 (z) ln(M gas ,5oo/M 5 oo) = 
-2.37 + 0.21 ln(M 5 no/2 x 10 14 M Q ) to derive the cluster mass 
and radius (rsoo) from our gas mass estimate. 

2.2.2. X-ray luminosity 

The X-ray luminosity is estimated by integrating the X-ray sur- 
face brightness. At 3<x significance, the surface brightness pro- 
files are detected out to at least 7*500 f° r a ll 62 clusters combin- 
ing XMM-Newton and ROSAT data (see Zhang et al. 2009). In 
practice, we estimate the total count rate from the background- 
subtracted, flat-fielded, point-source-subtracted, and point- 
spread-function (PSF) corrected surface brightness profile in 
the 0.7-2 keV band, and convert this to X-ray luminosity using 
the best-fit "mekal'0 model given by the spectra in XSPEC in 
the aperture covering all annuli defined in Sect. 3.2 in Zhang 
et al. (2009). We note that we do not study the temperature 
scaling relations for this sample here because of the inhomoge- 
neous range of projected distances used to measure the cluster 
temperature. 

We show the XMM-Newton-ROSAT vs. ROSAT-only mea- 
sured X-ray luminosity in the 0.1-2.4 keV band in Fig. IA.ll in 
Appendix^] The XMM-Newton-ROSAT to ROSAT-only mea- 
sured X-ray luminosity ratio is (92 + 2)% with (0.07 + 0.01) dex 
scatter. The faint point sources subtracted from the XMM- 
Newton data may account for a small fraction of the difference. 
A systematic difference in the flux calibration between ROSAT 
and XMM-Newton might play a major role in the 8% differ- 
ence in the X-ray luminosity (e.g., Snowden et al. 2002). A 
good fraction of the scatter may be introduced by the varying 
amounts of point sources and, especially, substructures that get 
excluded in the ROSAT and XMM-Newton analysis. 

In addition, there are some low-temperature systems (i.e., 
NGC507, Fornax, NGC1550, MKW4, NGC4636, NGC5044, 

2 The critical density is given by p c (z) = E 2 (z)3Hg(SnGy l , where 

eHz) = n m (i + z) 3 + n A + (i - a m - n A )(i + z f. 

3 For A0401, A0478, A1736, A2163, and A2244, which have no 
XMM-Newton data, we extracted the temperature from Vikhlinin et al. 
(2009a) assuming 0.3Z o iron abundance. 
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and A3581; kT < 2 keV) in the sample. To examine 
whether the line emission becomes important and boosts 
the X-ray luminosity for those systems, we show the 
iron abundance vs. temperature relation for the 62 clus- 
ters in Fig. IB. II in Appendix [B] The best fit is Z/Z = 
10 -(0.323±o.o6i) ( k r/keV) -(0.324±o.098) using the bisector method 

and Z/Z = 10-(0-325 ± 0.043) (fer/keV) -(0.320±0.068) using ±e QJ ._ 

thogonal method, respectively. This is consistent with the re- 
sults found in Balestra et al. (2007) but for clusters at higher 
redshifts (z > 0.3) and in a higher temperature range (3- 
15 keV), though their clusters show a steeper slope than that 
for our nearby clusters. The iron abundance vs. temperature 
correlation indicates that a flux-limited sample tends to include 
low-mass systems with high iron abundance, of which the X- 
ray luminosity is in part boosted by the line emission. This may 
modify the scaling relations at the low-mass end in terms of the 
mass dependence of the slope and the intrinsic scatter. 

Fabian et al. (1994) pointed out that some clusters are sig- 
nificantly above the best fit of the luminosity scaling rela- 
tion because of the presence of cool cores. This motivates the 
cluster core correction in deriving the X-ray luminosity (e.g., 
Markevitch 1998). We focus on the results using the X-ray lu- 
minosity corrected for the cluster core (hereafter L co ) by as- 
suming a constant value in the cluster core equal to the value 
at 0.2r 500 , S X (R < 0.2r 500 ) = Sx(0.2r 50 o) (Zhang et al. 2007). 
We note that this correction is only applied in determining the 
X-ray luminosity, not the gas mass. The bolometric luminos- 
ity corrected for the cluster core is listed in Table [1] and the 
bolometric luminosity within rgoo (hereafter L ln ) and in the 
[0.2- 1.0]r5oo radial range (hereafter L ex ) are listed in Table lC. ll 
To examine the scatter in the scaling relations caused by the 
presence of cool cores, we also compare the results using L co 
with those using L m (Appendix ICt and L ex (Appendix iDb. re- 
spectively. 

Since the soft band X-ray luminosity is widely used in stud- 
ies of the scaling relations, we calibrate the luminosity scal- 
ing relations using both the bolometric luminosity in the 0.01- 
100 keV band (Lboi) and the soft band luminosity in the 0.5- 
2 keV band (Lo.5-2keV, see also Appendix[E]i. 

2.3. Quantification of the cluster dynamical state 

2.3.1 . Offset between the X-ray flux-weighted center 
and BCG position 

The X-ray flux-weighted center of each cluster is listed in 
Cols. 2-3 of Table [T] which is determined based on XMM- 
Newton data as described in § 2.3 in Zhang et al. (2010). 
Our choice of the BCG position (Cols. 4-5 of Table [TJ is 
explained in § 12.11 The angular separation between the X- 
ray flux-weighted center and BCG position is converted into 
the physical separation at the cluster redshift, and is listed in 
TableQ] in units of both kpc (ofoffset) and r 500 (afoffsetMoo)- 

The offsets between the X-ray flux-weighted centers and 
BCG positions for the 62 clusters closely follow a log-normal 
Gaussian distribution (Fig. [2] left panel). The best fit of 
log 10 (£/ ff se t/r5oo) gives a mean value of -(1.93 + 0.06) and 
cr = (0.50 + 0.06). Forty-six clusters show < 0.037r 50 o off- 



sets, within lcr of the mean value. The remaining 16 clusters 
are sparsely spread over the range of [0.037 - ljrsoo- The best 
fit of log 10 (af o ff set /kpc) gives a mean value of (1.03 + 0.06) and 
cr = (0.55 + 0.06). Forty-seven clusters show < 38 kpc off- 
sets, within lcr of the mean value. The remaining 15 clusters 
are sparsely spread over the range of [38 - 1000] kpc. Thirteen 
of those 16 clusters with large offsets between the X-ray flux- 
weighted centers (see Table[TJ and BCG positions are disturbed 
clusters (see Tableland Sect. 12.3.31 . 

2.3.2. Central cooling time 

The central cooling time can be more accurately estimated 
from Chandra data because of its smaller PSF We thus use 
the central cooling time calculated at 0.004r5oo from Eq.(15) in 
Sect. 2.6 in Hudson et al. (2010) to divide the sample of the 62 
clusters into 26 cool-core clusters (i.e., "SCC" in Hudson et al. 
2010) and 36 non-cool-core clusters (i.e., "NCC" and "WCC" 
in Hudson et al. 2010) as listed in Table [2] Interestingly, we 
also found a correlation between the central cooling time and 
the offset between the cluster center and BCG position (Fig. [2] 
right panel). The best power-law fit to the relation between the 
offset and central cooling time is log 10 (^=^) = (-2.051 + 
0.058) + (0.907 ± 0.081) log 1Q ( Cool ^ time )Tnd log 10 (^) = 
(0.874 + 0.059) + (1.003 + 0.081) log 10 ( c °°'™s time ). 

2.3.3. X-ray morphology 

The combined MOS and pn images in the 0.7-2 keV band are 
shown in Appendix [F] According to their X-ray flux images, 
Vikhlinin et al. (2009a) divide the 62 clusters into 41 undis- 
turbed clusters and 21 disturbed clusters listed in Table|2] 

3. Results for the observational sample 

We investigate the three scaling relations between the lumi- 
nosity and velocity dispersion, luminosity and gas mass, and 
cluster radius and velocity dispersion, respectively, for the 62 
clusters in the HIFLUGCS. To examine possible systematic un- 
certainties due to the choice of the fitting method, we apply the 
BCES bisector and orthogonal methods. For all 62 clusters, the 
best power-law fits of all studied relations given by the bisector 
and orthogonal methods are consistent (Table [3]). We therefore 
focus on the best fits given by one of the two methods, i.e., the 
BCES bisector method, to illustrate the results. 

3.1. L- cr relation 

We summarize the best power-law fits of the L - cr relations 
using the X-ray bolometric luminosity (L^° ol ) and 0.5-2 keV 
luminosity (-£ u 5 _2i a .v)> respectively, in Table [3] In Fig. [3] we 
show the L c °. - cr relation of the 62 clusters. 

bol 

The slope for the 62 clusters, i.e., (4.02 + 0.33), agrees with 
the self-similar prediction (L ~ cr 4 ). The slopes for the undis- 
turbed clusters, disturbed clusters, cool-core clusters, and non- 
cool-core clusters are statistically indistinguishable. Ignoring 
their measurement uncertainties, both slopes for disturbed and 
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undisturbed clusters are steeper than for the combined sample. 
This is because most disturbed clusters are below the best-fit 
relation, most undisturbed clusters are above, and hardly any 
low-mass clusters are flagged as disturbed. The slope for the 
combined sample is thus influenced by a number of low-mass 
systems, which are all undisturbed clusters. The normalization 
for the undisturbed clusters is ~ 60% higher than for the dis- 
turbed ones. 

The intrinsic scatter (Table [3} of the undisturbed clusters 
and cool-core clusters is only slightly larger than that of the 
disturbed clusters and non-cool-core clusters. The clusters with 
more morphological substructure do not show larger scatter 
than those with less substructure. This indicates that the scat- 
ter driven by the presence of cool cores is comparable to that 
driven by substructure using L co . The increasingly large scatter 
toward the low-mass end is caused by the systems with < 45 
cluster members with spectroscopic redshifts. In Sect. [5] we 
will discuss the systematic uncertainties in the velocity disper- 
sion measurements due to the limited number of cluster mem- 
bers. 

The top-left panel of Fig.|4]shows the histogram of residu- 
als in logarithmic space from the Lboi.soo - cr relation. The best 
Gaussian fit gives 0.33^ ^ dex scatter, dominated by the intrin- 
sic scatter, i.e., (0.27 + 0.03) dex. We note that the histogram 
does not closely follow a symmetric Gaussian distribution, 
which may slightly underestimate the scatter. The top-right, 
bottom-left, and bottom-right panels of Fig.|4]show the residu- 
als as a function of the offset between the X-ray flux-weighted 
center and BCG position, luminosity fraction within 0.2r5oo, 
and central cooling time, respectively. There are very weak cor- 
relations caused mainly by the systems that have fewer than 
45 cluster members with spectroscopic redshifts, for which the 
measurement uncertainties in the velocity dispersion can be 
large and in part account for the scatter. 

As shown in Fig. [3] the normalization of the L — cr relation 
for the sample in Wu et al. (1999) is slightly higher than that 
of our sample. For L ln , the two samples are in better agreement 
(see Fig. IC.ll in Appendix ICb. Therefore, the core correction 
applied when deriving the X-ray luminosity for our sample ac- 
counts for the normalization difference between our sample and 
the sample of Wu et al. (1999) in Fig. [3] The different slopes 
between two samples may be due to their different selection 
functions as the sample in Wu et al. (1999) is not a flux-limited 
sample. 

The presence of cool cores is one of the main causes of the 
scatter in the L - cr relation as the scatter using L m for the non- 
cool-core clusters is ~ 20% smaller than that for the cool-core 
clusters (Appendix [Cj- When the X-ray luminosity corrected 
for the central region (< 0.2rsoo) is used, the intrinsic scatter 
is smaller by ~ 0.05 dex equaling (0.27 + 0.03) dex for the 
sample of 62 clusters. The intrinsic scatter in the L - cr rela- 
tion is similar using L c ° and L ex (Appendix iDt. The residuals 
of Lj". - cr are more strongly correlated with the luminosity 
fraction within 0.2r5oo than the residuals of both Lj^ - cr and 
- cr. The best-fit relation is A lgL = (0.52 + 0.10)°+ (0.51 + 
0.09) log 10 (Lboi,o.2r 500 /^bouoo) with its correlation coefficient of 
0.44. Correcting or excluding the central emission therefore ef- 
ficiently reduces the intrinsic scatter. 



3.2. L - M gas relation 

We summarize the best power-law fits of the L - M gas relations 
using the X-ray bolometric luminosity (IS) and 0.5-2 keV 
luminosity (io° g _2keV^ respectively, in Table [3] In Fig. [5] we 
present the L™j - M gas relation of the 62 clusters. The slope 
of the best-fit power-law for the 62 clusters is (1.29 + 0.05). 
The slopes for the undisturbed and disturbed clusters are statis- 
tically identical. The slope for the non-cool-core clusters, i.e., 
(1.42 + 0.05), is steeper than that for the cool-core clusters, i.e., 
(1.24 ± 0.06). The intrinsic scatter of those subsamples is com- 
parable. 

The top-left panel of Fig. [6] shows the histogram of resid- 
uals in logarithmic space from the best-fit I{j° 1500 - Afgas.soo 
relation for the 62 clusters using the BCES bisector method. 
The best Gaussian fit gives (0.07 + 0.01) dex scatter in loga- 
rithmic space, comparable to the intrinsic scatter. We note that 
the histogram does not closely follow a symmetric Gaussian 
distribution, which may slightly underestimate the scatter. The 
top-right, bottom-left, and bottom-right panels of Fig. [6] show 
the residuals as a function of the offset between the clus- 
ter X-ray flux-weighted center and BCG position, luminosity 
fraction within 0.2r5oo, and central cooling time, respectively. 
We do not observe as clearly evident correlations as for the 
^boi 5oo _ ^gas,500 relation in Fig. IC.ll in Appendix|C] This indi- 
cates that the cluster core correction may sufficiently suppress 
the scatter caused by the presence of cool cores. 

The intrinsic scatter in logarithmic space of the L - M gas 
relation using L co is similar to that using L ex (AppendixlDt. but 
is 0.05 dex lower than that using L m (AppendixICb. The cluster 
core correction in deriving the X-ray luminosity significantly 
reduces the intrinsic scatter in the L - M gas relation. 

We note that both quantities are derived from the X-ray 
surface brightness distribution in the soft band. If the gas 
is clumped, the emission measure can be overestimated by 
(«g)/(n e ) 2 , which results in overestimation of both the X-ray 
luminosity and the gas mass. Therefore, one possibly underes- 
timates the scatter in the relation. 

3.3. r 50 o - cr relation 

In Fig. [7] we present the relation between the velocity disper- 
sion and cluster radius, the latter being determined from the 
mass vs. gas-mass relation (see Sect. 12.2. U . In Table [3j we 
summarize the best power-law fits. The slopes for the undis- 
turbed clusters, disturbed clusters, cool-core clusters, and non- 
cool-core clusters are statistically indistinguishable. 

Surprisingly the intrinsic scatter in logarithmic space 
(Table [3]) for the cool-core clusters is about a factor of two 
larger than for the non-cool-core clusters. Since most undis- 
turbed clusters are cool-core clusters, the undisturbed clusters 
exhibit significantly larger intrinsic scatter than the disturbed 
clusters. This indicates that the presence of cool cores is the 
main driver of the scatter instead of the morphological sub- 
structure. We note that the scatter becomes increasingly large 
toward the low-mass end, which is coincidentally again caused 
by the systems that have fewer than 45 cluster members with 
spectroscopic redshifts. 
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The top-left panel of Fig. [8] shows the histogram of the 
residuals in logarithmic space from the rsoo - cr relation. 
The best Gaussian fit gives (0.063*qq19) scatter in logarithmic 
space, comparable to the intrinsic scatter. We note that the his- 
togram has a strong asymmetric shape, such that the Gaussian 
distribution slightly underestimates the scatter. The top-right, 
bottom-left, and bottom-right panels of Fig.[8]show the residu- 
als as a function of the offset between the X-ray flux-weighted 
cluster center and BCG position, the luminosity fraction within 
0.2r5oo, and the central cooling time, respectively. The residuals 
are not very tightly correlated with any of these three param- 
eters. However, 20 of the 26 cool-core clusters are above the 
best fit, and two thirds of the non-cool-core clusters are below 
the best fit. The undisturbed and disturbed clusters display ho- 
mogeneously distributed residuals. This confirms that the pres- 
ence of cool cores is the main cause of the intrinsic scatter in 
the rsoo - cr relation. 

4. Simulated vs. observational samples 

To understand the cluster physics behind the observed scaling 
relations, it is crucial to compare observational samples to rep- 
resentative samples in simulations. Our analysis of our sam- 
ple shows that the presence of cool cores is one of the main 
causes of the scatter. In addition, it has become increasingly 
clear that active galactic nuclei (AGN) play an important role 
in understanding the properties of clusters (e.g., McNamara & 
Nulsen 2007) and their scaling relations. We therefore investi- 
gated how well simulations can explain the observed results by 
comparing our observational measurements to those for a sam- 
ple of 2 1 clusters and groups simulated at a very high resolution 
both with and without AGN feedback (Puchwein et al. 2008). 
The AGN feedback model that was employed resolves some 
of the long-standing problems that hydro-dynamical simula- 
tions of galaxy clusters typically have, i.e., excessive overcool- 
ing within the densest cluster regions and too bright and too 
blue central galaxies. The AGN feedback model also brings the 
simulated X-ray luminosity - temperature scaling relation into 
excellent agreement with the observational one. We note that 
simulations with different cluster physics (e.g., Borgani et al. 
2004; Evrard et al. 2008) may give different predictions about 
the normalization, slope, and scatter of the scaling relations. 

4.1. A sample drawn from simulations 

Puchwein et al. (2008) carried out a set of high-resolution 
hydrodynamical re-simulations of clusters selected from the 
Millennium simulation with and without AGN feedback (for 
the AGN feedback model see Sijacki et al. 2007), and present 
the corresponding L - T relation. In Puchwein et al. (2010), 
they also show the properties of the stellar components and 
halo baryon fractions for the same sample. Their simulations 
have high enough resolution to accurately resolve galaxy pop- 
ulations down to the smallest galaxies that are expected to con- 
tribute significantly to the stellar mass budget. We select a sam- 
ple of 21 galaxy clusters from their simulations, whose gas 
masses span a similar range as our observational sample, i.e., 
(2.95 x 10 11 - 1.10 x 10 14 )M o for the case without AGN feed- 



back and (0.74 x 10 11 - 1.23 x 1O 14 )M for the case with AGN 
feedback. The gas masses for the same sample of 21 simulated 
clusters but with AGN feedback expand to a broader range due 
to the feedback. 

4.2. Analyzing the sample from simulations 

Both X-ray bolometric luminosity (also corrected for the clus- 
ter core) and velocity dispersion are derived in the same manner 
as for the observational sample of the 62 clusters except that 
we do not use the caustic method to identify interlopers since 
they are known in simulations. In determining the velocity dis- 
persion, we exclude sub-halos that do not contain any stellar 
components and should thus not be considered as galaxies. 

In the simulations, the virialized region of every cluster 
is known. We derive both the velocity dispersion using those 
galaxies within a projected radius of 1.2 Abell radii (hereafter 
crdiity), and the velocity dispersion using those galaxies not only 
within a projected radius of 1 .2 Abell radii, but also within the 
virialized region of the cluster (hereafter cr c ] ean fl In the cr c i ean 
case, nearby galaxies and interlopers are completely removed. 
Therefore, cr clean gives a reliable estimate of the velocity dis- 
persion. We note that even cr dirty contains only interlopers that 
are very close to the cluster center since they are all within the 
high resolution region of the cluster re-simulation, which typi- 
cally extends to five times the cluster virial radius or somewhat 
farther depending on its exact geometry. 

4.3. L™ - cr relation 

The L£°j - cr dirty relation is slightly shallower than the L£°j - 
cr c iean relation for the simulated sample (Fig. [3j. This may sug- 
gest that interlopers bias the velocity dispersion estimates to- 
ward lower values, which might be the case for the six clusters 
that have fewer than 45 cluster members with spectroscopic 
redshifts in the observational sample. We discuss this further in 
Sect. [5] The scatter in the L-cr relation of the simulated sample 
is larger for the run with AGN feedback. 

The shape and scatter of the L-cr relation of the sim- 
ulated sample with AGN feedback is comparable to that of 
the cool-core clusters in the observational sample, disregard- 
ing the six clusters with < 45 cluster members with spectro- 
scopic redshifts. AGN feedback suppresses cool cores, thus re- 
duces the X-ray luminosity in simulations. AGN feedback also 
significantly lowers halo gas mass fractions in low-mass sys- 
tems. We therefore find that the scatter becomes larger toward 
the low-mass end, i.e., systems with cr < 500 km s , which 
is also present in our observations. The observational sample, 
disregarding the six clusters with < 45 cluster members with 

spectroscopic redshifts, gives a best fit of log 10 ( WT^ L ^ ) = 



4 More precisely, we include only galaxies that are part of the 
friends-of-friends (FoF) group of the cluster in calculating cr cltm . 
The FoF group is computed using a linking length of 0.2 times the 
mean inter-particle distance, and extends to roughly ~ 1.7r2oo which 
corresponds to ~ 1.2 Abell radii for clusters with a virial mass of 
2.7 x 10 14 M Q /A. For more massive clusters, the FoF group is of course 
larger, while it is significantly smaller for poor groups. 
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(4.46 + 0.23) log 10 (j^rr) + (31.40 + 0.66), which closely fol- 
lows the simulated sample with AGN feedback. 

In the high-mass regime, the observational sample shows 
that non-cool-core clusters are the main driver of the scatter. 
Since most non-cool-core clusters are disturbed clusters, their 
substructures cause overestimations of the cr values (see Fig. 6 
and Fig. 9 in Biviano et al. 2006). Different fractions of sub- 
structures therefore translate into scatter in the L - cr relation, 
which is exactly what we find in the observational sample. 
The scatter in the L - cr relation of the simulated sample is 
smaller than that of the observational sample in the high-mass 
regime since the simulated sample does not predominantly con- 
tain mergers. In part, this difference might also be due to too 
few massive clusters in our simulations, i.e., six systems with 
0"ciean > 500 km s _1 . Forthcoming simulations of much larger 
cosmic volumes will be very useful in differentiating the scat- 
ter in the L - cr relation caused by measurement systematics 
from that attributable to cluster physics and achieving a clearer 
understanding of the cluster dynamics and gas physics in the 
high-mass regime. 

4.4. L - M gas relation 

The intrinsic scatter in the L-M gas relation is small for both the 
observational sample and the simulated sample. AGN feedback 
mainly tends to move clusters downward along the L - M gas 
relation rather than strongly changing the relation because re- 
moving gas from within rsoo also significantly reduces the X- 
ray luminosity. Hence, AGN feedback does not produce sig- 
nificant scatter in the L - M gas relation found in simulations. 
Nevertheless, the simulated sample with AGN feedback has a 
slightly steeper slope than that of the case without AGN feed- 
back. The difference between the simulated samples without 
and with AGN feedback is still small, and comparable to the in- 
trinsic scatter. There is a good agreement between the L - M g . ds 
relations from observations and simulations as shown in Fig. [5] 
Interestingly, the best fits of the disturbed clusters and undis- 
turbed clusters in the observational sample closely match the 
simulated samples with and without AGN feedback, respec- 
tively. This suggests that one of the main causes of the scatter 
could be AGN activities providing feedback in different phases, 
e.g., during a feedback cycle. 

4.5. r$oo - cr relation 

As shown in Fig. [7] the slope of the rsoo _ c relation for the 
observational sample is similar to that of the simulated sam- 
ple with AGN feedback, disregarding the clusters again that 
have fewer than 45 cluster members with spectroscopic red- 
shifts. For the simulated sample, the fraction of gas removed by 
AGN feedback becomes significant toward low-mass systems. 
As a consequence, their DM distributions expand slightly. Both 
the removal of gas and the expansion of the DM distribution re- 
sult in increasingly smaller cluster radii with decreasing mass 
in the simulated sample with AGN feedback compared to the 
simulated sample without AGN feedback. The simulated sam- 
ple with AGN feedback thus has a slightly steeper slope than 



that of the sample without AGN feedback. For the observa- 
tional sample, the subsample of the undisturbed clusters ex- 
hibit a steeper slope than that of the subsample of the disturbed 
clusters. This is consistent with the scenario that incorporates 
AGN activities in the undisturbed clusters (mostly cool-core 
clusters). 

There is some offset in the normalization between the sim- 
ulated sample and the observational sample, which is how- 
ever still within the scatter. X-ray masses are lower than the 
true masses in numerical simulations (e.g., Evrard 1990; Lewis 
et al. 2000; Rasia et al. 2006; Nagai et al. 2007; Piffaretti 
& Valdarnini 2008; Jeltema et al. 2008; Lau et al. 2009; 
Meneghetti et al. 2010). This may in part account for the off- 
set in the normalization, which relies on the total mass vs. gas 
mass calibration. The galaxy selection is complete for the sim- 
ulated cluster. However, we have no homogeneous photometry 
data to constrain the completeness for the clusters in the obser- 
vational sample. Differences in the selection of galaxies used 
to compute the cr may also in part cause this offset. 

5. Systematic errors in estimates of cr 

5.1. Galaxy selection by projected radial distance 

Most velocity dispersion profiles of galaxy clusters become flat 
beyond l/i _I Mpc which suggests that the measured velocity 
dispersion within a larger radius is more representative of the 
total kinetic energy of the cluster galaxies (e.g., Fadda et al. 
1996; den Hartog & Katgert 1996; Biviano & Girardi 2003; 
Boschin et al. 2010). 

We test how the radial selection of cluster members affects 
the velocity dispersion estimates of the observational sample 
as follows. We measure the velocity dispersion within [0.5, 
1.0, 1.5, 2.0, 2.5]xr5oo, and normalize it to the value mea- 
sured within 1.2 Abell radii (top-left panel of Fig. [§}. On av- 
erage, the velocity dispersion measured within small radii, i.e., 
[0.5,1.0]xr5oo, is ~ 10% larger than the one measured within 
larger radii. This is consistent with den Hartog & Katgert 
(1996) finding that more clusters with relatively large veloc- 
ity dispersion than small when measuring velocity dispersion 
close to the cluster center. We also note that the scatter in the 
measured velocity dispersion within small radii is ~ 3 times 
that measured within 2.5rsoo- 

In the top-right panel of Fig. [9] we show the normalized 
velocity dispersion measurements within rsoo as a function of 
the velocity dispersion measured within 1.2 Abell radii. For 
systems of velocity dispersion greater than 500 km s , there 
is on average less than 10% difference between the velocity 
dispersion measurements within rsoo and within 1 .2 Abell radii. 
The difference becomes larger for low-mass systems, and is up 
to ~ 30% on average for our sample. The scatter in the ratio of 
the velocity dispersion measurements within rsoo and to those 
within 1.2 Abell radii is almost independent of the absolute 
value of the velocity dispersion, at ~ 25%. 

In the bottom panels of Fig. [9] we also show the normalized 
velocity dispersion measurements within rsoo as a function of 
the offset between the X-ray flux-weighted cluster center and 
BCG position and the central cooling time. For clusters with a 
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smaller offset between the X-ray flux-weighted cluster center 
and BCG position or shorter central cooling time, the veloc- 
ity dispersion measurements within rsoo are significantly larger 
than the values measured within 1.2 Abell radii. 

We note that interlopers introduce uncertainties in the 
above tests, particularly for the six systems with < 45 cluster 
members with spectroscopic redshifts. We therefore also car- 
ried out tests using the simulated sample as shown in Figs.fTOl- 

m 

The trends of velocity dispersion decrease with increasing 
radius agree between the simulated sample and the observa- 
tional sample. In Fig. [10] the simulated sample shows that AGN 
feedback does not clearly affect the velocity dispersion esti- 
mates. However, interlopers increase both the amplitude and 
the scatter in the deviations of the <x estimates as a function of 
projected cluster-centric distance. The average deviation for the 
simulated sample without interlopers is similar to that of the 
observational sample. However, the scatter for the simulated 
sample with interlopers is comparable to that of the observa- 
tional sample. This indicates interlopers may affect the velocity 
dispersion estimates for a few but not the majority of systems 
in the observational sample. 

As shown in Fig. [10] the velocity dispersion within 
1 .2 Abell radii for two groups in the simulated sample is ~ 30% 
larger than that within smaller cluster-centric radii. One group 
is in a strongly clustered region with several group-size objects 
within 1.2 Abell radius in projection. In particular, one of the 
group-size objects has a similar mass to the group we analyzed. 
The other group is in the process of merging, which biases the 
cr estimate toward larger values (see also Biviano et al. 2006). 

In Fig. Q~T] the simulated sample without interlopers con- 
firms that systems with velocity dispersions greater than 
500 km s~' have < 10% difference between the velocity disper- 
sion measurements within rsoo and within 1.2 Abell radii. The 
results for the simulated sample also indicates that interlopers 
can boost the scatter in the bias for the low-mass systems of 
velocity dispersion < 500 km s -1 . For the low-mass systems, 
the uncertainties in the velocity dispersion estimates may in- 
deed be as large as 40%, disregarding the radial selection. We 
have to keep this in mind when we consider the six clusters 
with < 45 cluster members with spectroscopic redshifts in the 
observational sample. 

5.2. Mass selection 

For the observational sample, we collected cluster galaxy red- 
shifts from the literature. This may introduce a bias in the ve- 
locity dispersion estimates because brighter cluster galaxies 
may be more likely to have published redshifts than fainter 
ones. This bias becomes less significant when many cluster 
galaxies with spectroscopic redshifts are available. Since our 
observed sample of cluster galaxies is incomplete, we test how 
the mass selection of cluster members affects the velocity dis- 
persion estimates using the simulated sample, which is homo- 
geneous in terms of cluster galaxies. In Fig. [12] we display the 



velocity dispersion determined for a fraction of cluster mem- 
bers at the massive end for the 2 1 simulated cluster^ 

AGN feedback does not have an obvious effect on the ve- 
locity dispersion estimates. For the simulated sample, the local 
regression non-parametric fit illustrates that the velocity dis- 
persion estimate tends to be increasingly biased toward smaller 
values with decreasing fraction of cluster members. The bias 
on average is within a few per cent as long as more than 10% 
of the cluster members at the massive end are used. The scatter 
in the velocity dispersion also increases as a smaller fraction 
of cluster galaxies are used, and is within 10% when at least 
50% cluster members at the massive end are used. The scat- 
ter is slightly smaller when there are no interlopers. As shown 
in Figs. IG.ll4G.2l when we consider only 45 of the most mas- 
sive cluster members, the uncertainties in the velocity disper- 
sion estimates can be up to 40% for some low-mass systems 
(cr < 50010ns- 1 ). 

5.3. Interlopers 

Except for one system in the simulated sample, interlopers al- 
ways bias the measurements of the velocity dispersion toward 
smaller values (see also Biviano et al. 2006). A significant frac- 
tion of galaxies (up to ~ 50% of n ga i) within 1.2 Abell radii are 
not in the virialized region for poor systems. This is not the case 
for massive systems. As shown in Fig. [13] a caustic cannot effi- 
ciently exclude interlopers at larger radii, i.e., [1 -2.5]r5oo, and 
may significantly bias the measurements of the velocity disper- 
sion toward smaller values for poor systems. The cr c i ean is a far 
more robust indicator of the cluster mass than the cr dirty for poor 
systems. 

We note that cr dirty for the simulated sample only contains 
interlopers very close to the cluster. In the observations, there 
may be more distant interlopers for poor systems. As shown in 
Fig. 10 in Biviano et al. (2006), unrecognized interlopers that 
are outside the virial radius but dynamically linked to the host 
cluster and do not form major substructures, bias the cr estimate 
toward smaller values than cluster galaxies. 

6. Conclusions 

We have presented the L - cr, L - M eas , and r^m - cr relations 
for the 62 clusters in the HIFLUGCS, a purely X-ray flux- 
limited sample selected to minimize bias toward any cluster 
morphology. The systems in this sample span a broad range 
of morphological substructure, central cooling time, and off- 
set between the X-ray flux-weighted cluster center and BCG 
position, respectively. Owing to our representative, statistically 
large sample, with ~ 1 .3 Ms of clean X-ray XMM-Newton data 
and 13439 spectroscopically confirmed cluster members for 62 
clusters, we have been able to minimize our measurement un- 
certainties in both X-ray observables and velocity dispersion. 
Our main results are as follows: 



5 Note that only galaxies above the resolution limit of the simula- 
tions are included in Fig.fl2]and Figs, IG.fUG.2l The three figures are 
meant to illustrate trends instead of giving quantitative constraints. 
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- The luminosity vs. velocity dispersion relation agrees with 
the self-similar prediction. The presence of cool cores is 
one of the major contributors to the scatter in the L - cr re- 
lation. Correcting the central region in deriving the X-ray 
luminosity reduces the intrinsic scatter from 0.33 dex to 
0.27 dex. Even after correcting the X-ray luminosity for the 
cool core, the scatter caused by cool cores becomes increas- 
ingly large toward the low-mass end. The scatter caused 
by the non-cool-core clusters does not strongly depend on 
the mass range, but becomes dominant for massive sys- 
tems. The intrinsic scatter for the non-cool-core clusters, 
0.25 dex, is statistically indistinguishable from that of the 
cool-core clusters, 0.28 dex, after correcting the central re- 
gion when deriving the X-ray luminosity. 

- The presence of cool cores is also one of the major con- 
tributors to the scatter in the L - M gas relation. Using the 
X-ray luminosity corrected for the cool core, the disturbed 
clusters with significant X-ray substructures exhibit simi- 
lar scatter as the undisturbed clusters, partly because of the 
preponderance of cool-core clusters in the undisturbed sub- 
sample. 

- The shape of the L c ° - cr relation in simulations with AGN 
feedback matches the observational sample, specifically the 
cool-core clusters, disregarding the clusters that have fewer 
than 45 cluster members with spectroscopic redshifts. A 
common trend in both observations and simulations is that 
the scatter becomes larger toward the low-mass end, i.e., 
systems with cr < 500 km s -1 . The shape and intrinsic 
scatter in the L co - cr relation of the observational sam- 
ple closely matches that of the simulated sample for the 
low-mass clusters indicating that AGN feedback operates 
there. In the high-mass regime, the observational sample 
shows that non-cool-core clusters (their substructures) are 
the main driver of the scatter. The scatter in the L - cr re- 
lation at the high-mass end is larger than the scatter in the 
simulated sample. This may be in part because there are too 
few massive clusters and no significantly disturbed clusters 
in the simulated sample. 

- Interestingly, the best fits of the luminosity vs. gas mass re- 
lations for the disturbed clusters and undisturbed clusters 
in the observational sample closely match those of the sim- 
ulated samples with and without AGN feedback, respec- 
tively. This suggests that one of the main causes of the scat- 
ter could be AGN providing feedback in different phases, 
e.g., during a feedback cycle. 

- The rsoo - cr relation of the observational sample is simi- 
lar to that of the simulated sample, disregarding the clus- 
ters with < 45 cluster members with spectroscopic red- 
shifts. For the simulated sample, the fraction of gas re- 
moved by AGN feedback becomes significant toward low- 
mass systems, which makes their potential wells shallower. 
The slope for the simulated sample with AGN feedback is 
thus steeper than that for the sample without AGN feed- 
back. For the observational sample, the subsample of the 
undisturbed clusters exhibits a steeper r 5 oo - cr relation than 
that of the subsample of the disturbed clusters. This sug- 
gests that there is AGN activity in the undisturbed clusters, 
which are mostly cool-core clusters. 



- Both the selections of the aperture and mass limit of the 
cluster members and interlopers cause systematic uncer- 
tainties in estimating the velocity dispersion. For the ob- 
servational sample, the scatter in the velocity dispersion 
measured within small radii, i.e., [0.5,1.0]xr 5 oo is ~ 3 
times that measured within 2.5r5oo- The analysis of the 
simulated sample indicates that interlopers bias the veloc- 
ity dispersion estimates toward smaller values. The inter- 
lopers increase both the amplitude and the scatter in the 
bias, which is particularly significant for low-mass systems 
(cr < 500 km s _1 ). The scatter in the bias of the veloc- 
ity dispersion estimates increases as the fraction of clus- 
ter galaxies used decreases. The scatter is slightly smaller 
when there are no interlopers. 
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Table 1. Offset between the X-ray flux-weighted cluster center and BCG position, velocity dispersion, and X-ray bolometric 
luminosity. 
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Table 2. XMM-Newton observations and cluster properties. 
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"S" denotes cool-core clusters, and "W" and "N" denote non-cool-core clusters. 
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Table 3. Power-law fit, log 10 (F) = A + Blog 10 (X), to the scaling relations of the observational sample. 
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Fig. 1. Line-of-sight velocity vs. projected radius (left panels) and sky positions (right panels) of the selected galaxies in a rich 
cluster, i.e., Coma (top panels) and in a poor cluster, i.e., SI 101 (bottom panels). 
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Fig. 3. X-ray bolometric luminosity vs. velocity dispersion with luminosity corrected for the cluster core Our observational 

sample is shown in red (undisturbed) and blue (disturbed) colors, with filled triangles and open boxes denoting cool-core and 
non-cool-core clusters. The black circles highlight the six clusters with < 45 cluster galaxy redshifts in the determination of the 
velocity dispersion. The black, red, and blue lines are the best fits using the BCES bisector method for the whole observational 
sample, subsample of the undisturbed clusters, and subsample of the disturbed clusters, respectively. The simulated sample is 
shown in black (with AGN feedback) and gray (without AGN feedback) stars using cr^ny. Crosses show the corresponding cases 
using <r c iean> the velocity dispersion being based only on those galaxies within the virialized region of the cluster and within a 
projected radius of 1 .2 Abell radii for the simulated sample. It is worth noting that no redshift correction and cool-core correction 
is applied in Wu et al. (1999). 
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Fig. 4. Top-left: Histogram of residuals in logarithmic space from the best-fit Lj^ - cr relation for the 62 clusters using the BCES 
bisector method. Top-right: Residual vs. offset between the X-ray flux-weighted center and BCG position. Bottom-left: Residual 
vs. fraction of the X-ray luminosity within 0.2rsoo- Bottom- right: Residual vs. central cooling time. The colors and symbols have 
the same meaning as those in Fig. [3] 
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Fig. 5. X-ray bolometric luminosity vs. gas mass with luminosity corrected for the cluster core The colors, lines, and 

symbols have the same meaning as those in Fig. [3] 
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Fig. 6. Top-left: Histogram of residuals in logarithmic space from the best-fit L™ t - M gas relation for the 62 clusters using the 
BCES bisector method. Top-right: Residual vs. offset between the X-ray flux-weighted center and BCG position. Bottom-left: 
Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Bottom- right: Residual vs. central cooling time. The colors and 
symbols have the same meaning as those in Fig. [3] 
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Fig. 7. Gas-mass-determined cluster radius vs. velocity dispersion. The colors, lines, and symbols have the same meaning as 
those in Fig. [3] 
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Fig. 8. Top-left: Histogram of residuals in logarithmic space from the best-fit rsoo - cr relation for the 62 clusters using the BCES 
bisector method. Top-right: Residual vs. offset between the X-ray flux-weighted center and BCG position. Bottom-left: Residual 
vs. fraction of the X-ray luminosity within 0.2f5oo- Bottom- right: Residual vs. central cooling time. The colors and symbols have 
the same meaning as those in Fig. [3] 
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Fig. 9. Top-left: Velocity dispersion within a projected radius of R normalized by the velocity dispersion within 1 .2 Abell radii as a 
function of the projected radius. We do not include the values for the clusters having fewer than 10 members within the projected 
radius we are interested. Top-right: Normalized velocity dispersion within a projected radius of rsoo vs. velocity dispersion. 
Bottom-left: Normalized velocity dispersion within a projected radius of r^oo vs. offset between the X-ray flux-weighted center 
and BCG position. Bottom- right: Normalized velocity dispersion within a projected radius of 7*500 vs - central cooling time. The 
results are only based on the observational sample of the 62 clusters. The colors and symbols have the same meaning as those in 
Fig. [3] The curves are the local regression non-parametric fits. 
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Fig. 10. Velocity dispersion measured by the galaxies within a projected radius of R normalized by the velocity dispersion within 
1.2 Abell radii as a function of the projected radius for the simulated sample. The results are only based on the simulated sample 
of the 21 clusters. The colors and symbols have the same meaning as those in Fig. [3] The curves are the local regression non- 
parametric fits. We do not include the values for the clusters having fewer than 10 members within the projected radius we are 
interested. The black circles highlight the derived velocity dispersion with < 45 cluster members. 
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Fig. 11. Velocity dispersion measured by the galaxies within a projected radius of rsoo normalized by the velocity dispersion 
within 1.2 Abell radii as a function of velocity dispersion for the simulated sample. The results are only based on the simulated 
sample of the 21 clusters. The colors and symbols have the same meaning as those in Fig. [10] The curves are the local regression 
non-parametric fits. 
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Fig. 12. Velocity dispersion measured by the n most massive galaxies normalized by the velocity dispersion within 1 .2 Abell radii 
as a function of the fraction of galaxies for the simulated sample. The results are only based on the simulated sample of the 21 
clusters. The colors and symbols have the same meaning as those in Fig.[lO] The curves are the local regression non-parametric 
fits. 
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Fig. 13. Line-of-sight velocity vs. projected radius of the 30 brightest member galaxies (left panels) and all members within 
1.2 Abell radii (right panels), respectively, for a simulated cluster (without AGN feedback) having 40 cluster galaxies when 
excluding interlopers. The top panels correspond to cr dirty , and the bottom panels correspond to o- c i ean . 
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Appendix A: Luminosity cross-calibration 

To cross-calibrate the XMM-Newton-ROSAT with the ROSAT- 
only measured X-ray luminosity, we re-derived the X-ray lumi- 
nosity from ROSAT within r^oo given in Sect. 12.2. ll bv using the 
gas mass from the current work and the mass vs. gas mass rela- 
tion in Pratt et al. (2009). The same spectral model was used to 
derive the X-ray luminosity using both ROSAT data alone and 
a combination of XMM-Newton and ROSAT data. The compar- 
ison between the XMM-Newton-ROSAT and ROSAT-only mea- 
sured luminosity in the 0. 1-2.4 keV band is shown in Fig. IA.ll 
The XMM-Newton-ROSAT to ROSAT-only measured lu- 
minosity ratio is (92 + 2)%. The intrinsic scatter is (0.07 + 
0.01) dex. This was found for the REFLEX-DXL sample of 14 
massive galaxy clusters at z ~ 0.3 in Zhang et al. (2006) and the 
REXCESS sample of 31 nearby galaxy clusters in Pratt et al. 

(2009, -^-r = 1.15 x (-^)° 94 ). The difference between the 

v ' erg s 1 \ erg s 1 / ' 

XMM-Newton-ROSAT and ROSAT-only measured luminosity 
is well within the intrinsic scatter. 

Appendix B: Iron abundance vs. temperature 
Appendix C: Scaling relations using L m 

In Table IC.ll we present the X-ray bolometric luminosity 
within rsoo (Z^" j). We list the best fits to the corresponding scal- 
ing relations using the bolometric and 0.5-2 keV band lumi- 
nosity in Table lC.2l and show those plots using the bolometric 
luminosity in Figs. IC.ll4C.2l which helps us to understand the 
scatter driven by the presence of cool cores. 
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use this image to demonstrate the existence of morphologi- 
cal substructure in the cluster. Significant substructure features 
shown in the image are excised before we perform the spectral 
and surface brightness analysis. 

As addressed in Sect. 12.31 13 of the 16 clusters with large 
offsets between the X-ray flux-weighted centers (see Table Q} 
and BCG positions are disturbed clusters (see Table|2]i. We now 
comment on these 13 clusters. The BCGs in A0399 and A1736 
are slightly offset from the main X-ray emission. The ICM in 
A3376, A0754, A2256, and A3667 exhibits a comet-like tail, 
and their BCGs are at the opposite end from the X-ray centers 
probably because of their on-going dynamical activity. A3395s 
is the south component of a bi-cluster, and its BCG is at an 
X-ray weak bright peak. The ICM in A 1367 has multi-peaks, 
and the BCG is at the northwest X-ray peak, which is not the 
brightest one. The BCG in A2163 (A2255) is not a dominant 
BCG, which sits slightly east (west) of the X-ray center. This 
also applies but less significantly to some more clusters in the 
sample. The ICM in Coma, A3558, and A2065 shows some 
weakly disturbed features, and their BCGs are only 40-60 kpc 
away from the X-ray centers. A3158, A3391, and A0576 are 
relaxed clusters, and their BCGs are xL 40 kpc away from the 
X-ray centers. 

Appendix G: Systematic errors in estimates of a 



Appendix D: Scaling relations using L ex 

Since the luminosity derived in the [0.2-1] radial range 
is widely used to reduce the scatter caused by the presence of 
cool cores, we also present the X-ray bolometric luminosity in 
the [0.2-1] r 50 o radial range (L£*) in Table |CT] We also list 
the best fits to the corresponding scaling relations using the 
bolometric and 0.5-2 keV band luminosity derived in the [0.2- 
1] 7-500 radial range in Table ID. fl and show the plots using Lr\ 
inFigs. |P~TUD~2l 



Appendix E: Scaling relations using ^o°5-2kev 

We present the corresponding scaling relations using the 0.5- 
2 keV band luminosity corrected for the cluster central regions, 
Lg° 5 2keV , in Figs. IE.ll4E.2l The best fits are listed in Table[3] 



Appendix F: XMM-Newton images of the sample 

As the soft band is insensitive to the cluster temperature and 
has data of high signal-to-noise ratio, we use the MOS and 
pn combined image in the 0.7 - 2 keV band to illustrate the 
X-ray morphological substructure of each cluster (Figs. IF. 14 — 
IF.51 >. X-ray point-like sources are identified and subtracted. The 
holes, where the point-like sources were, are re-filled with the 
Chandra CIAO routine "dmfilth" using randomization based 
on the surface brightness distribution around the holes. We only 
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Table C.l. X-ray bolometric luminosity within /-500, L m , and in the [0.2 - ljrsoo annulus, L ex . 
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Note: The bolometric X-ray luminosity within the cluster core can be derived by L b oi,o.2r 500 = - £^ 
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Fig.A.l. XMM-Newton-ROSAT vs. ROSAT-only measured luminosity in the 0.1-2.4 keV band within rsoo. The dashed line 
denotes 1:1. With a fixed slope to 1 , the best-fit normalization of the XMM-Newton-ROSAT vs. ROSAT-only measured luminosity 
for the 62 clusters is 0.92 shown in solid line. The colors and symbols have the same meaning as those in Fig. [3] 
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Fig. B.l. Iron abundance vs. temperature for the 62 clusters. The black line denotes our best power-law fit using the bisector 
method. The dot-dashed line is the best fit in Balestra et al. (2007) for clusters at higher redshifts (z > 0.3) and in a higher 
temperature range (3-15 keV). The colors and symbols have the same meaning as those in Fig. [3] 
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Table C.2. Power-law fit, log 10 (F) = A + B log 10 (X), to the scaling relations for the observational sample using L" 
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Fig. C.l. Upper panel: X-ray bolometric luminosity vs. velocity dispersion with luminosity derived from all emission interior to 
r 5Q0 C^toi)- Lower left panel: Histogram of residuals in logarithmic space from the best-fit LP. - cr relation for the 62 clusters 
using the BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux- weighted center and BCG position. 
Lower 3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. central cooling 
time. The colors, lines, and symbols have the same meaning as those in Fig. [3] 





V 0ffset ( r 500) L bol,0.2.r 50 / L bol,500 Coolin 9 time ( G V r ) 

Fig. C.2. Upper panel: X-ray bolometric luminosity vs. gas mass with luminosity derived from all emission interior to rsoo C^i)' 
Lower left panel: Histogram of residuals in logarithmic space from the best-fit LS. - M gas relation for the 62 clusters using the 
BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux-weighted center and BCG position. Lower 
3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. central cooling time. The 
colors, lines, and symbols have the same meaning as those in Fig. [3] 
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Table D.l. Power-law fit, log 10 (F) = A + B log 10 (X), to the scaling relations for the observational sample using L e 
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Fig. D.l. Upper panel: X-ray bolometric luminosity vs. velocity dispersion with luminosity derived from emission in the [0.2 - 
1.0]r5oo aperture (IS). Lower left panel: Histogram of residuals in logarithmic space from the best-fit ZS - cr relation for the 
62 clusters using the BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux-weighted center and 
BCG position. Lower 3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. 
central cooling time. The colors, lines, and symbols have the same meaning as those in Fig. [3] 
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Fig. D.2. Upper panel: X-ray bolometric luminosity vs. gas mass with luminosity derived from emission in the [0.2 - l.OJrson 
aperture Lower left panel: Histogram of residuals in logarithmic space from the best-fit Lf^ - M gas relation for the 62 

clusters using the BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux-weighted center and BCG 
position. Lower 3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. central 
cooling time. The colors, lines, and symbols have the same meaning as those in Fig. [3] 




Fig. E.l. Upper panel: X-ray luminosity in the 0.5-2 keV band vs. velocity dispersion with luminosity corrected for the cluster 
core (£o°5_2kev)- Lower left panel: Histogram of residuals in logarithmic space from the best-fit L^° 52keV _ °~ relation f° r the 62 
clusters using the BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux-weighted center and BCG 
position. Lower 3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. central 
cooling time. The colors, lines, and symbols have the same meaning as those in Fig. [3] 
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Fig. E.2. Upper panel: X-ray luminosity in the 0.5-2 keV band vs. gas mass with luminosity corrected for the cluster core 
(Lg° 5 2keV ). Lower left panel: Histogram of residuals in logarithmic space from the best-fit £o° 5 _ 2 kev _ ^ gas re l at i° n f° r tne ^2 
clusters using the BCES bisector method. Lower 2nd panel: Residual vs. offset between the X-ray flux-weighted center and BCG 
position. Lower 3rd panel: Residual vs. fraction of the X-ray luminosity within 0.2r5oo- Lower right panel: Residual vs. central 
cooling time. The colors, lines, and symbols have the same meaning as those in Fig. [3] 




Fig. F.2. Combined MOS and pn images of the clusters in the 0.7-2 keV band, where point sources have been excised and refilled 
with values from neighboring pixels. 




Fig. F.3. Combined MOS and pn images of the clusters in the 0.7-2 ke V band, where point sources have been excised and refilled 
with values from neighboring pixels. 




Fig. F.4. Combined MOS and pn images of the clusters in the 0.7-2 keV band, where point sources have been excised and refilled 
with values from neighboring pixels. 




Fig. F.5. Combined MOS and pn images of the clusters in the 0.7-2 keV band, where point sources have been excised and refilled 
with values from neighboring pixels. 
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Fig. G.l. Velocity dispersion measured by the 45 most massive galaxies normalized by the velocity dispersion within 1.2 Abell 
radii as a function of the fraction of galaxies for the simulated sample. The results are only based on the simulated sample of the 
21 clusters. The colors and symbols have the same meaning as those in Fig.[lO] The curves are the local regression non-parametric 
fits. 
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Fig. G.2. Velocity dispersion measured by the 45 most massive galaxies normalized by the velocity dispersion within 1 .2 Abell 
radii as a function of velocity dispersion for the simulated sample. The results are only based on the simulated sample of the 21 
clusters. The colors and symbols have the same meaning as those in Fig. [10] The curves are the local regression non-parametric 
fits. 



